#!/usr/bin/python
# starcounts.py
'''Contains data about observed star counts from Cox 1999.'''

# Galactic latitudes for which data is available
bvals = [0, 5, 10, 20, 30, 60, 90]
# Magnitudes for which data is available
MvVals = [12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25]
# List of [lower bound, upper bound] pairs derived from MvVals
limitset = [[11.5, 12.5],[12.5, 13.5],[13.5, 14.5],[14.5, 15.5],[15.5, 16.5],[16.5, 17.5],[17.5, 18.5],[18.5, 19.5],[19.5, 20.5],[20.5, 21.5],[21.5, 22.5],[22.5, 23.5],[23.5, 24.5],[24.5, 25.5]]

# Data, log10(stars deg**-2). Outer index is magnitude, inner index is latitude.
N_bM = [[2.466, 2.434, 2.279, 2.011, 1.799, 1.450, 1.374],
        [2.922, 2.877, 2.671, 2.355, 2.138, 1.816, 1.728],
        [3.316, 3.262, 3.023, 2.678, 2.461, 2.149, 2.045],
        [3.650, 3.594, 3.344, 2.991, 2.772, 2.445, 2.326],
        [3.951, 3.897, 3.653, 3.301, 3.069, 2.702, 2.575],
        [4.243, 4.191, 3.955, 3.594, 3.341, 2.926, 2.799],
        [4.548, 4.494, 4.247, 3.861, 3.584, 3.133, 3.005],
        [4.868, 4.807, 4.530, 4.103, 3.804, 3.327, 3.191],
        [5.191, 5.120, 4.802, 4.326, 4.005, 3.503, 3.353],
        [5.479, 5.398, 5.039, 4.515, 4.173, 3.649, 3.488],
        [5.710, 5.623, 5.234, 4.673, 4.313, 3.770, 3.604],
        [5.909, 5.816, 5.400, 4.806, 4.430, 3.874, 3.706],
        [6.092, 5.991, 5.546, 4.916, 4.528, 3.970, 3.797],
        [6.258, 6.150, 5.669, 5.002, 4.606, 4.055, 3.875]]

# Allsky average for same magnitudes
allsky = [2.001, 2.396, 2.753, 3.075, 3.376, 3.663, 3.947, 4.230, 4.510, 4.761, 4.967, 5.145, 5.308, 5.453]
